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This work examines the autonomous navigation accuracy achievable for a lunar explo- 
ration trajectory from a translunar libration point lunar navigation relay satellite, aug- 
mented by signals from the Global Positioning System (GPS). We also provide a brief 
analysis comparing the libration point relay to lunar orbit relay architectures, and discuss 
some issues of GPS usage for cis-lunar trajectories. 


Nomenclature 

L\ Cis-lunar libration point (between Earth and Moon) 

L 2 Trans-lunar libration point (beyond Moon’s orbit) 

L 4 , Z /5 Leading and trailing equilateral libration points, respectively 
•£ 1 , £ 2 1575 MHz and 1227 MHz GPS carrier frequencies, respectively 


Introduction 

Farquhar 1 described several libration point navigation concepts that would appear to support NASA’s 
current exploration vision. One concept is a lunar relay satellite operating in the vicinity of Earth-Moon 
1/2, providing “Eaxth-to-lunax fax-side and long-range surface-to-surface navigation and communications 
capability.” Reference 1 lists several advantages of such a system in comparison to a lunar orbiting relay 
satellite constellation. Among these are one or two vs. many satellites for coverage, simplified acquisition 
and tracking due to very low relative motion, much longer contact times, and simpler antenna pointing. An 
obvious additional advantage of such a system is that uninterrupted links to Earth avoid performing critical 
maneuvers “in the blind.” Renault and Scheeres 2 have estimated that the stationkeeping cost for such a 
satellite, when considering the impact of navigation errors of 10 km and 10 mm/sec, would be on the order 
of 2-3 m/sec per month, assuming maneuvers approximately every two days. Folta and Vaughn 3 present a 
survey of stationkeeping costs for libration point orbits in the Earth-Moon system. 

Another concept Farquhar described is the use of Earth-Moon L\ for lunar orbit rendezvous, rather 
than low lunar orbit as was done for Apollo. Farquhar claims this technique requires only slightly higher 
fuel cost than low lunar orbit rendezvous for short-stay equatorial landings. More recently, Condon and 
Wilson 4 have estimated that libration point rendezvous missions have significant advantages over lunar 
orbit rendezvous missions when global access, long stay times, and anytime aborts are driving requirements 
for human exploration missions. However, these trades, and the trade vs. lunar surface rendezvous, are 
significantly complicated by as yet undetermined abort requirements. 

Farquhar also described an interplanetary transportation system that would use libration points as ter- 
minals for an interplanetary shuttle. This approach would offer increased operational flexibility in terms of 
launch windows, rendezvous, aborts, etc. in comparison to elliptical orbit transfers. More recently, other 
works 5,6 have shown that patching together unstable trajectories departing Earth-Moon libration points 
with stable trajectories approaching planetary libration points may also offer lower overall fuel costs than 
elliptical orbit transfers. The lunar navigation infrastructure should evolve to support such concepts. 

“Aerospace Engineer, Flight Dynamics Analysis Branch, Code 595, AIAA Senior Member, 
t Aerospace Engineer, Flight Dynamics Analysis Branch, Code 595, AIAA Member. 


1 of 15 


American Institute of Aeronautics and Astronautics 



Another concept Farqhuar described was a deep space relay at Earth- Moon L 4 and/or L 5 that would serve 
as a high data rate optical navigation and communications relay satellite. The advantages in comparison 
to a geosynchronous relay are minimal Earth occultation, distance from large noise sources on Earth, easier 
pointing due to smaller relative velocity, and a large baseline for interferometry if both La and L s are used. 
Such a relay could initially support lunar missions as well. 

Barton et al. 7 studied the use of the Global Positioning System (GPS) for navigation enroute between the 
earth and the moon. Assuming modest modifications that would improve GPS receiver sensitivity by approx- 
imately 10 dB and a high-gain directional receiver antenna, they showed that GPS signals viewed over the 
earth’s limb would support post-translunar injection (TLI) navigation out to about half the lunar distance. 
They also showed GPS navigation could support a mid-course trim burn for at least several hours after TLI, 
but if the trim bum was more than 8 hours after TLI, there was not enough GPS information to estimate the 
post-bum state. This level of GPS coverage might support the L\ lunar rendezvous scenario, especially if 
augmented by additional signals from NASA’s Tracking and Data Relay Satellite System (TDRSS), or from 
navigation assets in the vicinity of the moon. More significant improvements in GPS receiver sensitivity 
could likely be achieved by telemetering the GPS navigation message via a distinct communications link to 
lunar navigation receivers, avoiding the need for the users to decode the GPS ephemerides. Such a system 
might support limited GPS navigation all the way to the Earth-side lunar surface. 

In this work, we make a preliminary assessment of a lunar navigation infrastructure based on the concepts 
described by Farquhax and Barton et al. First we compare and contrast these concepts to various alternative 
navigation architectures for cis-lunar operations. Next, we describe calculations that establish the feasibility 
of using GPS for cis-lunar transfers, and possibly in support of L\ lunar rendezvous scenarios. Finally, 
we illustrate some representative autonomous navigation accuracy results for a cis-lunar transfer scenario 
that uses GPS pseudoranges in combination with one-way Doppler measurements from an Earth-Moon L 2 
Orbiter (EML20). 


Architecture Analysis 

In this section we describe and compare lunar-orbiting constellations and EML2 orbiting constellations. 
We did not perform an exhaustive analysis to determine optimal constellation coverage. Instead, we selected 
what appear to be representative circular and elliptical single-coverage constellations from the literature as 
a starting point for our lunar orbiting constellations, and a representative class of large amplitude orbits 
about Z/2‘ We also did not try to identify ‘frozen” lunar orbits, since these rely on properties of the lunar 
gravity field, and significant portions of the far-side field are currently poorly quantified due to the lack of 
lunar far-side tracking. 

Methodology 

In developing these concepts, we limited the trade space with the following guidelines: no outages during 
critical maneuvers and operations; navigation solutions must be computable onboard all user vehicles; navi- 
gation assets should be disposable at the end of mission; surface user latitude should be at least 75 degrees; 
minimum surface user elevation angle must be at least 10 degrees; a relatively small number of users must be 
supported in the lunar vicinity at any given time; and finally, instantaneous positioning is not required. We 
also assumed that an earth return vehicle remaining at L\ during a surface stay can serve as an Earth-side 
communications relay if needed to avoid long light-time delays that would result from using near-Earth relays 
for low-latitude Earth-side operations. Clearly, an additional relay satellite orbiting L\ could also serve this 
function if needed. We further assumed that all concepts would make maximal use of the Earth-orbiting 
GPS and TDRSS constellations. We did not explicitly study ground-based tracking from Earth since it may 
not support high-latitude surface operations and cannot support trans-lunar (far-side) operations, although 
it could clearly support cis-lunar space operations and near-side surface operations. 

Figure 1 depicts our reference mission, which consists of departure from a low Earth parking orbit at 
28.5 degree inclination, cis-lunar coast with lunar swingby to L\ with at least one TLI trim burn, parking 
and stage separation at L\, and descent to a high-latitude landing site. This reference mission has been 
selected not only as representative of a possibly typical exploration scenario, but also because it contains 
what we believe to be the events and mission phases that will most stringently drive navigation requirements: 
the TLI bum and post- TLI trim burns, L\ station-keeping, separation and rendezvous in the L\ vicinity, 
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Figure 1. Representative cis-Iunar trajectory viewed from above the moon’s orbit around the earth, in coor- 
dinates rotating with the moon’s orbit; this work primarily examines the segment labeled “1.” 


descent and landing to a high-latitude site, and surface operations at a high-latitude site. In this work, we 
focus primarily on the segment labeled u l n in the figure, which is the initial transit from Earth to Moon. 

Lunar Orbit Constellations 

As a starting point, we examined the possibility of establishing a four-satellite constellation that places three 
of the four spacecraft in high altitude, eccentric orbits. Constellations of this type, discussed by Draim 8 
for use around the earth, provide single-fold global coverage with a minimum number of satellites. The 
constellation Reference 8 proposes calls for three spacecraft with eccentricities of 0.13 and inclinations of 
33 degrees uniformly distributed in right ascension so as to yield maximum coverage of the globe. With the 
periapses for such orbits in the southern hemisphere, this is an ideal means of providing continuous coverage 
of the land masses of the earth, which axe concentrated in the northern hemisphere. For complete global 
coverage a fourth spacecraft is set into an equatorial orbit with a period half that of the other three. The 
semi-major axes of the three inclined spacecraft are determined geometrically by the need to maximize the 
coverage. If we scale this by the radius of the central body, we find that the 106,440 km semi-major axis 
called for at the earth translates into a semi-major axis of 29,000 km at the Moon. The problem here is that 
the coverage geometry scales with distance while the gravitational attraction of the central body scales with 
mass, which is a function of volume. Where the earth with its relatively large mass is the central body, this 
is not a problem since third body perturbations from the moon and sun have little effect on the long term 
stability of the constellation. When this geometry is scaled down for the moon however, we find that the low 
gravity of the Moon combined with the relatively large third body perturbations from the earth prevent the 
orbits needed for a constellation of this type about the moon from being stable. Further we find that if we 
were to reduce the semi-major axes of the constellation to address this stability issue, gaps in the coverage 
open up. This eliminated this particular elliptical orbit constellation from further consideration for global 
coverage at the moon. 

Next, we looked at how six satellites in a constellation of a type discussed for Earth coverage by Lang 
and Meyer 9 could be applied to the problem of continuous global lunar coverage. Reference 9 claims that 
this constellation provides continuous global coverage by at least one satellite at lower altitudes for a given 
minimum elevation constraint than any other six satellite constellation. This constellation, which Figure 2 
illustrates, employs four spacecraft in circular orbits, all at 58.9 degree inclinations, distributed in right 
ascension at 90 degree intervals. Two spacecraft in opposing circular equatorial orbits fill the coverage gaps 
left by the first four. The result is a constellation with a global coverage pattern that repeats at a semi- 
major-axis-dependent frequency. When the minimum elevation for a user anywhere on the lunar surface is 
10 degrees, a Moon-centered constellation with a semi-major axis of 6214 km and a period of 12.2 hours 
results. Such low, circular orbits about the moon turn out to be stable in the long term but can be costly 
to establish. For computing the establishment cost, we assumed the six spacecraft were sent to the moon 
with two launches, each carrying three spacecraft. In this scenario, each triad would establish a lunar orbit 
before separation from the carrier spacecraft. 

Since the trade space of possible Moon-centered elliptical constellations is so expansive, we developed a 
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Circular (Lang-Meyer) Configuration 



Figure 2. Circular lunar orbit constellation, based on the work of Lang and Meyer. 


Elliptical (Modified Lang-Meyer) Configuration 



Figure 3. Elliptical lunar orbit constellation, based on extension of the Lang and Meyer orbits. 
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representative case. In this case, which Figure 3 shows, we employed the orientations of the Lang and Meyer 
constellation using eccentric rather than circular orbits. We set periapsis to the minimum radius that results 
in a greater than 10 degree minimum elevation for a user on the lunar surface. We set an outer boundary 
for apoapsis at 10,000 km; we observed that higher orbits tended to be less stable, but we did not perform 
an exhaustive analysis. The semi-major axis of the subsequent orbits turned out to be 8125 km with a 
period of 18.2 hours and an eccentricity of 0.23. Establishing the eccentric orbits may be accomplished by 
approaching the lunar surface at 100 km altitude and executing a lunar orbit insertion (LOI) burn that sets 
the apoapsis of the desired orbit at 9993.75 km. Once the spacecraft moves to apoapsis, a bum would be 
executed raising periapsis to its final value of 6256.25 km. We found that the cost of establishing such stable 
eccentric orbits in the Lang and Meyer orientations were marginally more expensive than the unstable high 
altitude circular orbits, but considerably less expensive than the highly stable low altitude circular orbits. 
Inclination changes done at apoapsis are much less expensive for the eccentric constellation than for the stable 
low altitude constellation. In this sense, this eccentric version of the Lang and Meyer constellation may be 
used as representative of the expense of establishing a useful, Moon-centered global coverage capability. Here 
again, we computed the establishment cost assuming that two TLI vehicles are used to get two triads of 
spacecraft to the moon, which upon LOI separate to form the constellation. 

Earth-Moon L 2 Orbiter Constellation 

Next, we sought to define a trajectory for one or 
more Earth-Moon L 2 Orbiters that would provide 
coverage for high-latitude sites on the moon. Fig- 
ure 4 depicts one possible configuration, which is 
a compromise between minimizing station-keeping 
costs and maximizing the out-of-plane amplitude 
of the trajectory around L 2 so as to provide high- 
latitude coverage on the moon. This configuration 
consists of a pair of large Lissajous orbits that have 
opposite inclinations in the plane normal to the 
moon’s velocity. Although a single satellite in one 
of the orbits will provide continuous coverage for 
nearly the entire lunar far-side hemisphere, in order 
to get continuous coverage above 10 degree eleva- 
tion at both lunar poles requires two satellites in 
each orbit. Figure 5 shows the elevation angles to 
the constellation of EML2 Orbiters that a station located directly at either pole of a perfectly spherical 
Moon would see, ignoring the physical libration of the lunar pole (which is only a few thousandths of a 
degree). Table 1 lists the phases of the satellites in the two L 2 orbits. Note that these orbits have periods of 
about 15.5 days. Front side lunar coverage would be provided by near-Earth space networks such as TDRSS, 
Earth ground networks, or possibly by an Earth return vehicle stationkeeping near L\. For estimating the 
establishment cost, we assumed the four spacecraft would be sent to the moon with two launches, each 
carrying two spacecraft. 

Architecture Comparison 

Table 2 compares the architecture concepts we finally considered. Definitions of the table’s notation are as 
follows: n is the number of satellites required for global coverage, AF e is the establishment cost (including 
both Earth departure and final orbit injection), AF^ is the cost to dispose of the satellite, and L s is the 
maximum space loss, computed assuming maximum range from a polar site at zero elevation, and a frequency 
of 1.575 GHz. In deriving estimates for the AF, we did not consider navigation or maneuver execution errors, 
or finite burn losses. To estimate the disposal AF, we chose the lower cost between impacting the moon or 
escaping from the lunar vicinity (we expect that zero disposal cost will be associated with the L 2 concepts 
since by merely ending station-keeping, they will escape the lunar vicinity) . We did not draw any conclusions 
regarding stationkeeping costs, since this will depend very strongly on as yet undetermined requirements. 
Preliminary investigations indicate however that the lunar orbits we examined are stable over at least a 
one year timeframe. The libration orbits will obviously require stationkeeping maneuvers, likely as often as 


Table 1. Phases of the EML2 orbiters, which are de- 
scribed with respect to the maximum elevation above 
the lunar orbit plane of a virtual satellite at the initial 
simulation epoch of 17 April 2008, 0.0 UT. For this cov- 
erage pattern, satellites 1 and 3 are in the blue “upper” 
1/2 orbit, and satellites 2 and 4 are in the red “lower” 
orbit, which Figure 4 depicts. 


Satellite 

Phase [days] 

1 

1.94 

2 

-1.94 

3 

-5.81 

4 

5.81 
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Figure 4. Earth-Moon L 2 Orbiters’ trajectories. The direction of the moon’s motion is toward the viewer, 
and the earth is to the right of the scene. The amplitudes of the orbits in the direction normal to the moon’s 
orbit plane are approximately 25,000 km, and they extend approximately 50,000 km on either side of the 
Earth-Moon-L 2 line. The amplitude of the motion toward and away from the moon covers about 25,000 km. 
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Figure 5. Plan view of L 2 Orbiters’ trajectories along with elevations from stations at lunar poles, assuming 
alignment of Moon’s axis with Moon’s orbit’s angular momentum. For this coverage pattern, satellites 1 and 3 
are in the blue “upper” L 2 orbit, and satellites 2 and 4 are in the red “lower” orbit. The green dashed lines 
in the X — Z plot indicate 10 deg. elevations from the lunar poles. The discontinuities in the trajectories are 
locations where bi-weekly stationkeeping maneuvers would need to be performed. 
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every revolution (Le. on the order of every two weeks). Based on the work of Reference 3, we estimate a 
stationkeeping cost for these types of orbits at around 300 meters per second per year. In any case, whatever 
AV advantage the EML20 concept may have over lunar orbiters must be offset against a maximum space 
loss penalty of approximately 20 dB. However, the maximum space loss is less than 10 dB worse than for 
a geosynchronous Earth orbit (GEO) satellite, and the minimum space loss for the EML20 concept is only 
about 2 dB worse than for a GEO satellite. As a result of this brief and incomplete survey, we believe 
that a case can be made for lunar orbiting constellations in support of long-term human and robotic lunar 
exploration, but we believe that the EML2 Orbiter concept provides a great deal of capability with minimal 
investment as a starting point for a comprehensive lunar and planetary navigation and communications 
infrastructure beyond near Earth orbit. 

Table 2. Comparison of representative architecture options described in the text; units of AV are km/sec; 
units of space loss are dB. Note that AV e includes departure from low Earth orbit, and both AV columns 
represent rough estimates of the total cost for the entire constellation, exclusive of trajectory knowledge and 
control errors. 


Concept 

n 

A^e 

A V d 

L s 

Circular lunar orbits 

6 

13.7 

2.2 

-172 

Elliptical lunar orbits 

6 

11.4 

2.3 

-176 

Large lissajous at L 2 

4 

8.7 

0 

-194 


GPS Signal Availability Analysis 

In this section, we evaluate the GPS signal availability over a four day time span covering the Earth-Moon 
transit segment “1” in Figure 1. 

Methodology 

As a baseline case, we considered a GPS receiver with an acquisition/tracking sensitivity of 25 dB-Hz 
and a 10 dB receiving antenna. The 25 dB-Hz threshold corresponds to a minimum ambient power of 
approximately —184 dB Watts (dBW), or received power at the output of the antenna of —174 dBW. This 
performance is consistent with a state of the art GPS receiver implementing modem weak signal tracking 
techniques. 10 We also examined cases that assumed acquisition and tracking thresholds of 18 dB-Hz and 
11 dB-Hz. The former case is probably at or just beyond the limit of current state of the art weak-signal 
GPS technology to decode the GPS navigation message (broadcast ephemeris). The latter case is probably 
at or just beyond the limit of current state of the art weak-signal GPS technology to acquire and track the 
GPS signals at all. In this case, the receiver would require an alternative source of ephemeris information 
for the GPS satellites, which limits its autonomy. 

The modeled GPS link budget includes realistic antenna gain patterns for the GPS satellites and receiver, 
free space propagation losses, and antenna/receiver noise models to evaluate signal levels reaching the receiver 
and corresponding carrier to noise spectral density. Geometric obstruction of GPS signals with respect to 
the earth is evaluated based on a constant 6378.14 km Earth radius plus a 50 km atmosphere mask. A 
single frequency i\ coarse/acquistion code GPS receiver was assumed, and GPS transmitted power levels 
were based on the Block II/HA GPS satellites, transmitting at typical power levels, which are approximately 
1.5 dB above minimum specified levels. The GPS constellation was simulated based on a nominal 27 satellites, 
six plane configuration. 

Results 

Figure 6 illustrates representative GPS signal availability for the three cases described above. Figure 7 
shows the corresponding ambient power and carrier to noise ratio (C/Nq) levels for signals. Considering the 
baseline 25 dB-Hz case, there is nearly continuous GPS coverage with frequently two or more satellites in 
view during the first half day after the TLI burn, out to about 18 Earth radii {Re) from the earth. This 
provides important observations in support of a possible correction burn that would be performed in this 
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timeframe. Peak GPS signal levels axe approximately 38 dB-Hz, and the spacecraft is already at a distance 
from the earth above which any GPS side lobe signals (greater than 25 degrees off transmitting GPS satellite 
nadir) extend above 25 dB-Hz. Prom this point onward there will never be more than three satellites visible 
simultaneously, and typically no more than two. 

Figures 8 and 9 provide a closer look at the GPS signals available in the vicinity of the Earth-Moon L\ 
point, at a range of approximately 50 R E and 2.4 days (58 hours) elapsed time. For the 25 dB-Hz case 
(Figure 8) passes that are interrupted because the GPS satellite sets behind the earth are short, about 10-15 
minutes on either side of the earth’s limb. For GPS satellites in which the earth does not obstruct the line 
of sight, the passes are typically 50-60 minutes in duration. Peak signal levels do not exceed 31 dB-Hz, and 
there is usually only a single GPS satellite visible at one time. The longest outage (no satellites available) was 
approximately 90 minutes. Figure 9 shows the data for an 18 dB-Hz sensitivity. The increased sensitivity 
means side lobe signals are available at higher altitudes (although not out to L\), which ultimately reduces 
the duration of outages. In the vicinity of the moon, peak GPS signals levels are approximately 29 dB- 
Hz, but overall GPS signal availability is very similar to availability near L\. As a general approximation, 
path losses are approximately 15-16 dB larger than GEO in the region from the L\ point to the lunar 
surface. Referring back to Figure 6, note that the 11 dB-Hz case appears to provide nearly continuous GPS 
coverage all the way out to the lunar surface. These results may be overly optimistic as they are based 
on tracking GPS side lobe signals out to 60 degrees off the GPS satellite nadir, and the actual side lobe 
power may deviate significantly from the modeled pattern beyond 30-40 degrees off nadir. The improved 
GPS availability indicated in the 18 and 11 dB-Hz cases could nevertheless be realized using an even more 
directional receiving antenna. The modeled 10 dB antenna has a 3 dB beamwidth of about 20 degrees (half 
angle), but in the vicinity of L 1? all of the GPS signals are received within a 5 degree cone (half angle) 
centered around the earth. 

A number of factors were examined to assess the sensitivity of the GPS availability to some of the 
assumptions used above. The simulation was repeated using transmitted power levels and antenna gain 
patterns for the modernized Block HR-M GPS vehicles, and for the new civil “L2C” signal on the GPS 
^2 carrier. GPS coverage was not significantly different for either of these cases. The new £2 signal may 
actually have some advantages for weak signal tracking however, because it has a signal component that 
is not modulated with data allowing longer coherent integration times. 11 Additionally, the GPS coverage 
is dependent on the geometry of the GPS constellation with respect to the outbound trajectory of the 
receiver. A number of different outbound trajectories with different epochs were examined, and for the 
25 dB-Hz sensitivity, GPS satellites from only two of the six orbital planes were typically visible. In the GPS 
constellation, each plane may contain 4-6 satellites, so the peak number of satellites visible and maximum 
outages vary depending on which planes of GPS satellites are observed. 

Navigation Analysis 

In this section, we describe our methodology and results for obtaining representative estimates of the 
autonomous navigation accuracy for a user vehicle on segment “1” of the reference trajectory Figure 1 
depicts. 

Methodology 

We performed the navigation analysis using the GPS Enhanced Onboard Navigation System (GEONS). 
GEONS is an extended Kalman filter 12 with physically realistic process noise models 13,14 and a factored 
covariance matrix. 15 Its measurement model may include standard and singly-differenced GPS 16 and GPS 
Wide Area Augmentation System (WAAS) data, intersatellite crosslink data, one-way forward-link Doppler 
from ground stations or TDRSS, and line-of-sight and angular measurements to celestial objects 17 and other 
spacecraft. The GEONS state vector may consist of absolute or relative position and velocity vectors in 
the mean equator and equinox of J2000.0 coordinate system, 18 empirical drag and solar radiation pressure 
coefficient corrections, and measurement biases for all data types. GEONS uses the Joint Gravity Model-2 
(JGM-2), 19 with coefficients up to degree and order 30, and models solar, lunar, and planetary point mass 
gravity using either the Van Flandern 20 or Moshier 21 analytical methods. GEONS has an analytic repre- 
sentation of the Harris-Priester 22 atmospheric density model, and a simple solar radiation pressure (SRP) 
model, similar to that described in Vallado, 23 that assumes a cylindrical Earth shadow. The filter integrates 
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Figure 8. GPS signals above 25 dB— Hz near the Earth- Moon L\ point. 
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Figure 9. GPS signals above 18 dB— Hz near the Earth-Moon L\ point. 
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the equations of motion with either 4th- or 8th-order Runge-Kutta methods, and measured accelerations may 
be incorporated. The state transition matrix is propagated with a semi-analytic formulation that includes 
Earth point mass and J 2 , atmospheric drag, and solar radiation pressure partial derivatives. 

For this analysis, we simulated measurements with a truth ephemeris that we generated from initial condi- 
tions just after the cis-lunar injection departure maneuver from low Earth orbit. The truth model parameters 
were as follows: 30x30 JGM2 Earth gravity; point-mass Moon gravity; Moshier analytical ephemerides for 
the Moon, Sun, Venus, Mars, Jupiter, and Saturn; 4.51 square meters and 100 kilograms spacecraft area and 
mass; spacecraft drag coefficient of 2.0; spacecraft SRP coefficient of 1.1886, with 1371 Watts per square 
meter assumed as the solar flux; and Allan variance parameters for the receiver clock typical of an oven- 
stabilized crystal oscillator (ho = 2.4 x 10“ 22 sec, and h _2 = 8.0 x 10 -28 sec -1 ). To numerically integrate 
the equations of motion, the truth model used an 8th-order fixed-step Runge-Kutta integrator with 3rd-body 
forces and rotation matrices evaluated at each function call, with a time step of 60 sec. For the simulated 
GPS pseudorange measurements, we added white Gaussian noise with standard deviation of 2 meters for 
signals stronger that 40 dB-Hz, and 10 meters for weaker signals. We simulated GPS transmitter clock 
and ephemeris errors using the ‘Tear 4” 24 model with 2 meter standard deviation, and used a curve fit of 
an averaged electron density profile for ionosphere errors. To the 1575 MHz Doppler data from the EML2 
Orbiter, we added white Gaussian noise with standard deviation of 20 mHz. 

We used a somewhat simpler set of model parameters in the navigation filter: 8x8 JGM2 Earth gravity; 
zero drag; and Allan variance parameters for the receiver clock typical of an temperature-stabilized crystal 
oscillator (ho — 3 x 10 -14 sec, and h_ 2 = 8.4 x 10 -20 sec -1 ); all other parameters were the same as 
the truth model, except that we added zero-mean Gaussian random position, velocity, clock bias, and clock 
drift initial condition errors and a random walk SRP bias to the filter. The standard deviations of the 
initial errors were, for position: 31.6 meters in the radial and crosstrack directions, and 100 meters in-track; 
for velocity: 10 cm/sec in all directions; for clock bias and drift: 100 meters and 10 cm/sec; and for SRP 
bias: 0.01%. We assumed correlation coefficients of —0.95 between initial radial position and along-track 
velocity, and between initial along-track position and radial velocity. We computed a Gaussian random 
initial state perturbation for the filter using a Cholesky decomposition of the initial covariance. To account 
for dynamic model simplification, we also added random walk acceleration process noise, with root power 
spectral densities (m/sec 3 / 2 ) of 1 x 10 -5 radial and 1 x 10 -4 cross-track. For the SRP bias, the root noise 
density was 1 x 10 -5 sec -1 / 2 . The filter assumed a 10 meter standard deviation for the GPS pseudorange 
measurement noise, and a 600 mHz standard deviation for the EML2 Orbiter Doppler noise. The filter 
processed all available measurements once per minute, but it ignored GPS signals that passed less than 
1000 km above the earth’s limb. Filter tuning and monte carlo analysis of the filter performance were 
beyond the scope of our study, so one should view our results as representative of the kinds of accuracies 
that might be achieved, rather than as a definitive evaluation of the achievable performance. 

Results 

Figures 10 - 13 plot the results of the filtering investigation. In Figure 10, the upper subplot shows the GPS 
pseudorange residuals, surrounded by the three-sigma formal error envelope of the residuals derived from 
mapping the filter’s covariance. The middle subplot shows the one-way Doppler residuals from the EML2 
orbiter and its formal error envelope. Both of these subplots also show the filter’s residual ratios along 
the bottom edge, with a scale on the right. These figures indicate that the filter’s residuals are generally 
zero-mean and uncorrelated, that GPS pseudorange residuals are generally below a few tens of meters, 
and that EML20 Doppler residuals are generally below a few tens of milliHertz, with occasionally large 
Doppler residuals that exceed the filter’s six-sigma edit ratio test. We believe these large residuals are due to 
truncation differences between our quadruple-precision data simulator and the double-precision wordlength 
used in GEONS. The bottom subplot indicates the total number of contacts to both GPS and the EML2 
orbiter. Since the EML2 orbiter is always visible, this plot shows that after the first half-day, only zero to 
three three GPS satellites are visible. Note that the GPS ephemeris we used in the navigation analysis was 
based on an actual current broadcast ephemeris set, rather than the nominal ephemeris set used for the 
visibility analysis above. As a result, one may note differences in the visibility at any given time, although 
the overall visibility statistics do not significantly differ. 

Figures 11 and 12 show collections of position- and velocity-related error statistics. In both these figures, 
the first three subplots show the errors projected onto local coordinates that are aligned with the radial 
and cross-track directions. The “out-of-plane crosstrack” direction is co-linear with the two-body angular 
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Figure 10. Measurement information processed by the filter. 


momentum of the receiver, and the “in-plane crosstrack” direction completes a right-handed system with the 
other two directions. The fourth subplots show the position and velocity magnitude errors, and Figure 11 
also plots the semi-major axis error. Both these figures also plot a formal error envelope derived from the 
filter’s covariance matrix. Since these are single-case results, the value of the actual errors in predicting 
real-world performance is somewhat reduced. Instead, one may view the formal error envelope as having the 
value of a linear covariance analysis, that is, as guideline of the range of accuracy that might be achievable. 
In this sense, the accuracies of better than 1 km and 5 cm/sec axe quite promising in comparison with 
the trans-lunar post-maneuver results on the order of kilometers and dozens of centimeters per second that 
Beckman and Concha 25 achieved during Lunar Prospector using two-way Doppler from Deep Space Network 
34- and 26-meter tracking sites. It should be noted that the real-world results of Reference 25 were achieved 
in the face of all the constraints posed by an operational environment, including most notably imperfectly 
performed maneuvers that this study did not address. 

Figure 13 shows a collection of the bias error estimates of the filter. The first two subplots show clock 
bias and drift, and the lower subplot shows the correction to the SRP coefficient. These plots show that, 
despite the assumption of an OCXO, the clock errors axe comparable in size to the position and velocity 
errors. Although it is difficult to discern from the lower subplot, we did verify that the SRP correction is 
observable; however we did not find the filter to be very sensitive to even much larger unknown SRP biases. 
For this reason, we chose a filter tuning that allows the filter to consider the uncertainty in this parameter 
in its gain computation, but does not actually result in any significant updates to the parameter. 

We also examined cases with improved GPS receiver sensitivity, using both 18 dB-Hz and 11 dB-Hz 
acquisition and tracking thresholds. Although these cases significantly arrested the decline in GPS coverage 
as distance increased from the earth, we found that none of the GPS-only scenerios had consistently stable 
results beyond The primary difficulty with these cases is that the clock states become highly correlated 
with the radial error states, as a result of a combination of poor geometry and poor visibility. 
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Figure 11. Transient (left axis) and “steady-state” (right axis) single-case position and semi-major axis (SMA) 
errors of the filter. 
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Figure 12. Transient (left axis) and “steady-state” (right axis) single-case velocity errors of the filter. 
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Figure 13. Transient (left axis) and “steady-state” (right axis) single-case bias errors of the filter. 


Summary and Conclusions 

In this work, we have exa m ined the feasibility of an onboard navigation system for cis-lunar exploration, 
based on one-way Doppler from an Earth-Moon L<i orbiter and GPS pseudorange measurements. Our study 
indicates that accuracies of better than 1 km and 5 cm/sec may be feasible, which are quite promising in 
comparison with the trans-lunar post-maneuver results on the order of kilometers and dozens of centimeters 
per second achieved during Lunar Prospector using two-way Doppler from Deep Space Network 34- and 26- 
meter tracking sites. Since these are single-case results, the value of the actual errors in predicting real-world 
performance is somewhat reduced. Instead, one may view these results as a guideline for the range of accuracy 
that might be achievable. Although we believe that a case can be made for lunar orbiting constellations in 
support of long-term human and robotic lunar exploration, we believe that the EML2 Orbiter provides a 
great deal of capability with minimal investment as a starting point for a comprehensive lunar and planetary 
navigation and communications infrastructure beyond near Earth orbit. 
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